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Abstract 

In a trapped Bose-Einstein condensate, subject to the action of an alternating external field, 
coherent topological modes can be resonantly excited. Depending on the amplitude of the external 
field and detuning parameter, there are two principally different regimes of motion, with mode 
locking and without it. The change of the dynamic regime corresponds to a dynamic phase transi- 
tion. This transition can be characterized by an effective order parameter defined as the difference 
between fractional mode populations averaged over the temporal period of oscillations. The behav- 
ior of this order parameter, as a function of detuning, pumping amplitude, and atomic interactions 
is carefully analyzed. A special attention is payed to numerical calculations for the realistic case 
of a quadrupole exciting field and the system parameters accessible in current experiments. 
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I. INTRODUCTION 



Bose-Einstein condensation (BEC) in weakly interacting dilute Bose gases is a topic of 
high current interest and intensive investigations, both experimental and theoretical, as can 
be lnfarred form th e b oo k Q and review articles fl fl S, fl fl B, Q. At Suffi c,en tly low 
temperature, the weakly interacting gas is almost completely Bose-condensed and is well 
described by the Gross-Pitaevskii equation [lj. Then the system is in its ground state. One 
could pose a question whether nonground-state Bose-Einstein condensates could be realized? 

The possibility of creating nonground-state Bose-Einstein condensates was advanced in 
Ref. 9J. As is clear, such a condensate has to be nonequilibrium, since in an equilibrium 
system at zero temeprature, the state with the lowest energy would be always preferable. 
The most important point is that the atomic system would possess a discrete spectrum, or at 
least, would have finite gaps between different parts of its spectrum, such as occur in optical 
lattices. A straightforward way of obtaining a purely discrete spectrum is by confining 
atoms in a trapping potential. Then, in order to transfer atoms from their ground BEC 
state to another nonground level, the system has to be subject to an oscillatory external 
field. The field oscillation is to be close to the frequency of the transition between the 
ground and the chosen excited state of the trap. The states of the trapped coherent atomic 
cloud are given by the set of solutions to the stationary Gross-Pitaevskii equation, which 
are termed topological coherent modes. Due to the nonlinearity of the system, caused by 
atomic interactions, the spectrum is not equidistant and allows for the selection of the desired 
levels by means of the resonant modulation of an external field. It has been mathematically 
accurately proved ^, Q that this resonantly excited condensate can be well approximated 



by an effective two-level system. The analytical proof has been confirmed 



ll| by direct 



numerical simulations for the time-dependent Gross-Pitaevskii e quat ion. The topological 
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coherent modes have also been studied theoretically in Refs. 

generation of a nonground-state BEC has recently been observed in optical lattices 1 

The resonant system with nonground-state BEC, generated by an applied modulating 
field, possesses several unusual properties, as is reviewed in Ref. [19|. One of the interesting 
features is the existence of two essentially different regimes of motion. In one of them, 
the fractional mode populations are locked in the regions close to their initial conditions, 
while in another regime, the mode populations oscillate in the whole diapason between 
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zero and one. The transition between these two dynamical regimes represents a kind of a 



dynamic phase transition. It has been shown 
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22j that this dynamical transition is 



equivalent to a phase transition in an effective averaged system, accompanied by critical 



phenomena around the phase transition. According to Landau and Lifshitz [23|, an order 
parameter is a quantity defined in such a way that it takes non-zero values, positive or 
negative, in the unsymmetrical phase and is zero in the symmetrical phase. So, similarly 
to phase transitions in equilibrium systems, it is possible to define an order parameter 
characterizing the dynamic phase transition. The role of such an order parameter for the 
system with resonantly generated topological modes is played by the difference between the 
time-averaged mode populations for the ground and the excited states. 

The aim of the present work is twofold. First, we analyze the order-parameter behavior for 
a wide range of different values of the detuning, pumping, and interaction strengths. Second, 
we accomplish a numerical investigation of the order-parameter features for a setup typical 
for current experiments, so that to choose realistic conditions for experimentally realizing 
the discussed generation of the coherent topological modes and observing the dynamic phase 
transition between the mode-locked and unlocked regimes. 
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re main equations 



The structure of the paper is as follows. In Section [TTl we recal 

nn 

describing the resonant excitation of the coherent topological modes [9|, 121[, characterizing 
the coupling between the ground and excited states of BEC. In Section IHIt we present 
the results of our investigation for the behavior of the order parameter under the varying 
amplitude of the excitation field, detuning from resonance, and for different strengths of 
atomic interactions. In section IIVI we explore the case of a realistic alternating field that 
could be realized in experiment, for which we consider an oscillatory magnetic quadrupole 
field. We summarize our results in the concluding section . 

II. MAIN EQUATIONS AND DEFINITION OF MODEL 



The mathematical model, rather accurately describing a trapped BEC of a dilute Bose 
gas, is the well-known Gross-Pitaevskii equation [2], which is equivalent to the nonlinear 
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Schrodinger equation 

dt ' 2m 



<9$ h 2 
= ifi— , = V 2 + f/ 4rap (r) + A s |$| 2 , (1) 



where U trap (r) is the trapping potential, v4 s = 4ti(N — l)H 2 a s /m , with m being atomic 
mass; a s , the s-wave scattering length; and N, the number of atoms. The stationary wave- 
function $ can be written in the form 

$ n = $ n (r,t) = ^(r)e-^, (2) 

where ip n (r) is a solution to the stationary Gross-Pitaevskii equation 

H[ip n ]<p n = E n ip n . (3) 

Since the Hamiltonian in equation ([3]) is non-linear, the spectrum E n is not equally 
spaced Because of this, it is possible to couple the ground BEC state in the confining 
potential to another mode (labeled p) by applying an oscillatory field [SI], e.g., given by 

V p (r,t) = V(r)cos(wt), (4) 

where uo is the angular frequency characterizing the time variation of the field and V (r) is 
a spatially dependent amplitude. If this frequency is chosen to be close to the transition 
frequency between the two modes, then, as has been shown [si, Q, the total wave-function 
$ can be written as the sum of two modes 

<£ (r, t) = c (f) (r) e-^/ h + c p (t) ^ (r) e~^ h , (5) 

where the population rij of each mode is given by 

nj = \ Cj {t)\ 2 . (6) 

This implies that if we wish to know the temporal evolution of these mode populations, 
we need to find out how the coefficients Cj evolve in time. Substituting Eq. (jHJ) into the 
Schrodinger equation 

i <9$ (r t) 

H + V(r,t)]$(v,t)=ih?^V, (7) 



we obtain the rate equations for the coefficients Cq and c 



p; 



^ = -ia 0p n p co - ^(3e lAuJt c p , (8a) 



do i 

-± = -ta p0 n c p - -^e-^co. (8b) 
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An accurate derivation of these rate equations can be done ay , Q by employing the 
Krylov-Bogolubov averaging technique [24J. In Eq. (jHJ), we have the transition amplitude, 
defined as 

Oik = y J Wi (r) I 2 (%* (r) | 2 - \ Vj (r) | 2 ) dv, (9) 
the coupling amplitude caused by the external field, defined as 

P = \J ^{v)V{v)ip p {v)dv, (10) 

and the detuning Au = uj — cj p> o, which is the difference between the field frequency and the 
transition frequency, given by huj p $ = E p — E . Note that amplitude (Q is unsymmetric, 
i.e., not necessarily q:q p = a p o. It comes from the non-linearity introduced by interaction. 
The interaction will have different values depending on the occupied state and that may be 
the cause of the obtained unsymmetry in cijk- 

Our physical system is similar to a condensate in a double-well potential. Here, the 
coupling field characterized by the parameter (3 has the role of tunneling and there is the same 



nonlinear term associated with collisions. Instead of using the two- mode model [25J, we prefer 
to use the Gross-Pitaevskii formalism which predict the form of the spatial distributions of 
wave function for both, ground and excited, condensates. 

The system of equations (jHJ) can be solved numerically for each given set of parameters. 
In order to simplify the description, the following dimensionless parameters are introduced: 

a = — , b= , 5= , (11) 

a p o a p o a p o 



with the time scaling 



t' = a p0 t. (12) 



This system of equations allows us to find the temporal behavior of the mode amplitudes 
characterizing the fractional mode populations (jSJ) m Fig- CEJ we present the temporal 

evolution of the mode populations for a resonant field (5 = 0), the interaction amplitude 
a — 1, and different values of the pumping amplitude b. We assume that, initially, all atoms 
are in the ground-state BEC, i. e., co(0) = 1 and c p (0) = 0. 

We see that there are two different types of behavior of the population dynamics. For 
b < 0.5, the ground state population is always larger than that of the excited state. While 
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FIG. 1: Time evolution of the population fraction of the ground state (solid line) and the excited 
state (dashed line) in the presence of a resonant field, 5 = 0, and a = 1 for (a) b = 0.4; (b) b = 0.5; 
(c) b = 0.5001; (d) b = 0.6. 
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for b > 0.5, the time averages of both population modes are the same, which means that 
it is possible to transfer all atoms from the ground to the excited state. The same effect 
can be obtained for other values of a as we can see in Fig. [2J with a = 0.1 and a detuning 
5 = 0.45. 

As is seen, it is easier to detect a trapped BEC in an excited state, when the dynamics 
correspond to the mode-unlocked regime shown in the cases (c) and (d) of Figs. [HandEJ since 
in these cases the system stays longer in the excited state. In order to deeper understand 
how the excitation depends on the parameters of the system, it is useful to analyze the 
behavior of the order parameter, as defined in Refs. 
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III. ORDER PARAMETER FOR DYNAMIC TRANSITION 



The order parameter, characterizing the dynamic phase transition for the case of the 



resonant BEC, with the generated topological modes, can be defined [2l| . as the difference 
rj between the time-averaged fractional mode populations, 

t] = n -n p , (13) 

where rTj is the time average of over the whole oscillation period. This order parameter 
clearly shows us what kind of dynamic regime occurs. Observing the character of the mode 
oscillations, we note that in cases (a) and (b) of Fig. [Tj and Fig. El 77 > 0, and in the (c) 
and (d) cases, r\ = 0. 

We have calculated the behavior of rj as a function of b, for < a < 2 with different 
amplitudes of the resonant field. These results are shown in Fig. [31 from where we see that 
there is an abrupt change of the order parameter for a > 0.8, while for a < 0.8, 77 goes 
smoothly to zero. Also, for 0.8 < a < 1, rj becomes negative for some values of b. When 
the order parameter is negative, this does not mean that the fraction of the excited mode 
is always larger than the ground state population. It simply shows that the BEC on the 
average spends more time in the excited state, which would favor the detection of condensed 
atoms in that state. 

To better illustrate the peculiarities of the population evolution, we show in Fig. H] the 
mode populations no and n p , for a = 0.8, for different values of b. In all these cases we set 
the resonant external field (5 = 0). As presented in Fig. [3] and Fig. [U when a = 0.8 and 
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FIG. 2: Time evolution of the population fraction in the ground state (solid line) and the excited 
state (dashed line) in the presence of a field with the detuning 5 = 0.45 and a = 0.1 for (a) b = 0.2; 
(b) b = 0.2749; (c) b = 0.2750; (d) b = 0.4. 
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FIG. 3: Order parameter rj as a functions of b for different values of a with a resonant field 
excitation. 

b = 0.57, we have 7/ ~ —0.174. Then n p never reaches one, i. e., we do not get a completely 
excited BEC. For conditions presented in Fig. H](d), the excited state population is dominant 
for the major part of time. 

The behavior of the order parameter when varying the detuning for a fixed interaction 
strength a = 0.1 is presented in Fig. [51 which shows rj as a function of b for different values 
of the detuning. This figure, together with Fig. [3j characterizes the overall behavior of the 
order parameter under varying the system parameters a, b and S. 

The question that is yet left open is how and under what particular conditions the pro- 



posal of obtaining the dynamical phase transition for a trapped BEC [2l| could be realized 
in laboratory. To this end, we need to study the behavior of the order parameter rj for 
a realistic setup of the applied field geometry, field amplitude, and detuning, so that to 
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FIG. 4: Fractions of mode populations in the ground state (solid line) and the excited state (dashed 
line) in the presence of a resonant field, 5 = 0, and a = 0.8 for (a) b = 0.55; (b) b = 0.56; (c) 
b = 0.565; (d) b = 0.57. 
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FIG. 5: Order parameter n as a function of b for different values of 5 with a = 0.1. 

identify the critical field that changes the dynamic regime of the mode excitation. The 
specifications for the dimensionless parameters a, b, and 5 for a typical laboratory setup is, 
therefore, necessary. In the following section, we calculate these parameters for a magnetic 
quadrupole field, which corresponds to the experiment in progress in our group, employing 
a 87 Kb condensate. 

IV. MAGNETIC QUADRUPOLE AS AN EXCITATION FIELD 



The amplitudes a 0p , a p0} and (3 cannot be obtained as exact analytic expressions, since 
there are no exact analytic solutions for the Gross-Pitaevskii equation. We can, however, 



obtain rather accurate approximate solutions using the Optimized Perturbation Theory 



26|, 
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whose detailed description can be found, e.g., in Ref. [3|. In this case, approximate expres- 
sions for the ground state and three excited-state wave functions, for a harmonic trap 



are given by 



U trap = ±(u?r 2 + u,y), (14) 



/ 2 \ !/ 4 

/MOOO^OOOJ e -(u 000 r 2 +v 000 zl)/2l^ ^ 

Wo = «oio fe) ^ e W^^, (16) 



-'000 



-*001 



A..2 ,,2 \ V 4 „ 

4 %01%)1 \ ^ e -(«ooir 2 +^ooi^)/2i^ ^ 



«100_^10o\ ( 1 _ u r\ e -{u 100 r*+v 100 z*)/2ll 



)10 ° " [ ^JT V 11 



(18) 



where Z r = ^yH/m uj r is the oscillator length, and u nm k and t> nm k are the variational param- 
eters, defined by minimizing the corresponding energy of each mode. 

The alternating magnetic field in a magnetic trap can be created as an oscillatory 
quadrupole magnetic field of the form 

B(r, t) = (Axx + Ayy - 2Azz) cosut, (19) 

which corresponds to the field formed by a pair of coils operating in the anti-Helmholtz 
configuration. Then the potential V(r) in Eq. (4) is given by 



V(r) = -jl ■ B = g F m F fi B A Vr 2 + 4z 2 , (20) 

where ft is the atomic magnetic moment, gp is the gyromagnetic factor, mp is the magnetic 
quantum number, and \xb is the Bohr's magneton. It is easy to see that the sole topological 
mode, from the lowest modes given above, which could be coupled to the ground-state 
mode by the potential ( 1201) . is the radial dipole mode 0ioo- For this mode, we have made 
calculations of the order parameter rj as a function of the field-amplitude parameter A. 
We have chosen this parameter A because it is easily controllable in laboratory, since it is 



proportional to the current in the coils constructed to produce the desired field 

12 



27]. 



We set all parameters of the system according to the experiment in progress, where 87 Rb 
is employed, having mass mo = 0.144 x 10 -21 g and the scattering length a s ~ 6nm. The 
magnetic trap has u r = 2tt x 120 Hz and u z = 2tt x 24 Hz. The number of atoms is taken 
as N = 10 4 , which corresponds to a realistic experimental number. Figure M shows the 
obtained results for different detunings. In this case, the transition frequency between the 
ground state and the radial dipole mode is cuioo.o = 27r x 190 Hz. This figure clearly shows 
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FIG. 6: Order parameter rj for the exciting magnetic quadrupole field as a function of the field 
parameter A for different values of the field detuning. 

that there exist critical points in the behavior of rj, which depend on the detuning. In Fig. 
|6l when, for instance, Au = —2n x 2 Hz, the critical point happens for A = 0.096 G/cm. 
Three different situations are seen in Fig. [6j For Auo = —2n x 13 Hz, we have the situation 
similar to that for a = 0.1 and 5 = 0.40, shown in Fig. [51 when rj < 0. In this case, the 
critical point occurs when A ~ 0.119 G/cm. We also have the case when the critical point 
occurs for rj > 0, which is analogous to the behavior obtained for a = 0.1 and 5 > 0.45 and 
shown in Fig. [5j Here, for Auj = 2n x 13 Hz, the critical point happens for A ~ 0.068 G/cm. 
In the case of a negative detuning, Auo = —2tt x 36 Hz, there is no abrupt change in the 
behavior of rj, similarly to the case shown in Fig. [5], where a resonant field is considered. 

V. CONCLUSION 
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The dynamics of the resonant excitation of topological coherent modes of a Bose-Einstein 
condensate in a harmonic trap is analyzed. The behavior of the effective order parameter 
is investigated as a function of the external field amplitude for a variety of different system 
characteristics, such as the detuning from the resonance and interaction amplitudes. An 
important observation is that the detuning can compensate the influence of atomic interac- 
tions, producing in the case of a weak interaction (a = 0.1), the behavior similar to that 
obtained for a strong interaction (a = 1) and a purely resonant field. The pumping field, in 
the setup of the quadrupolar geometry was employed to study the feasibility of observing 
the dynamic phase transition for typical experimental conditions. The frequency variation 
around the resonance leads to the different behavior of the order parameter, depending on 
whether the detuning is red- or blue- shifted. The obtained results can be used for the ex- 
perimental realization of the dynamic phase transition in the process of resonant generation 
of topological coherent modes. 

Finally, we shall mention that several calculation using Bogoliubov theory performed by 
Dziarmaga and Sacha 28j] shown that the effect of quantum and thermal depletion could be 
a serious limitation on quantum coherence of atomic BEC. Second question is the necessity 
of an efficient atomic population transference. If both, quantum and thermal depletion, 
affect the process, we cannot guarantee that the atoms which populates the excited state 
are a Bose-Einstein condensate (a coherent state). Following Ruostekoski and Walls |29| . 
the effects of decoherence due to noncondensed atoms on BECs shows that purity decays 
fast. Hence, one must have the interaction parameter much smaller than the decoherence 
time scale associated with the nonlinear parameter. Because our analysis shown that the 
macroscopic population of the excited state is associated with small values of b, we believe 
the effects of decoherence will be small. 
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